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Abstract 

In  this  paper  we  develop  a  new  analytical  methodology  for  the  evaluation  of  the  outage 
probability  of  cooperative  decode-and-forward  (DF)  automatic-repeat-request  (ARQ)  relaying 
under  packet-rate  fading  (fast  fading)  channels,  where  the  channels  remain  fixed  within  each 
ARQ  transmission  round,  but  change  independently  from  one  round  to  another.  In  particular, 
(i)  we  derive  a  closed-form  asymptotically  tight  (as  SNR  — >  oo)  approximation  of  the  outage 
probability,  (ii)  we  show  that  the  diversity  order  of  the  DF  cooperative  ARQ  relay  scheme  is 
equal  to  2L  —  1,  where  L  is  the  maximum  number  of  ARQ  retransmissions,  and  (iii)  we  develop 
the  optimum  power  allocation  for  the  DF  cooperative  ARQ  relay  scheme.  The  closed-form  ex¬ 
pression  clearly  shows  that  the  achieved  diversity  is  partially  due  to  the  DF  cooperative  relaying 
and  partially  due  to  the  fast  fading  nature  of  the  channels  (temporal  diversity).  With  respect 
to  power  allocation,  it  turns  out  that  the  proposed  optimum  allocation  scheme  depends  only  on 
the  link  quality  of  the  channels  related  to  the  relay,  and  compared  to  the  equal  power  allocation 
scheme  it  leads  to  SNR  performance  gains  of  more  than  1  dB.  Numerical  and  simulation  studies 
illustrate  the  theoretical  developments. 

Keywords:  Automatic-repeat-request  (ARQ)  protocol,  cooperative  decode-and-forward 
(DF)  relaying,  outage  probability,  wireless  networks. 


1  Introduction 


Conventional  wireless  networks  involve  point-to-point  communication  links  and  for  that  reason  do 
not  guarantee  reliable  transmissions  over  severe  fading  channels.  On  the  other  hand,  cooperative 
wireless  networks  exhibit  increased  network  reliability  due  to  the  fact  that  information  can  be 
delivered  with  the  cooperation  of  other  users  in  networks  [1]-[12].  In  particular,  in  cooperative 
systems  each  user  utilizes  other  cooperative  users  to  create  a  virtual  antenna  array  and  exploit 
spatial  diversity  that  minimizes  the  effects  of  fading  and  improves  overall  system  performance. 
Cooperative  communications,  also  known  as  relay  channels,  was  first  introduced  in  [7]  in  which 
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a  three-way  channel  was  analyzed  based  on  the  capacity  region.  In  [8]-[10],  relay  channels  have 
been  analyzed  from  an  information-theoretic  point  of  view.  With  respect  to  practical/realistic 
cooperative  communication  protocols  for  wireless  networks,  past  literature  includes,  but  is  not 
limited  to,  the  work  in  [l]-[6],  [11,  12]  and  the  references  therein. 

Automatic-repeat-request  (ARQ)  protocols  for  wireless  communications  have  been  studied  ex¬ 
tensively  in  the  past  and  proved  themselves  as  efficient  control  mechanisms  for  reliable  data  packet 
transmissions  at  the  data  link  layer  [13]-[19].  The  basic  idea  of  ARQ  protocols  is  that  a  receiver 
requests  retransmission  when  a  packet  is  not  correctly  received.  Recently,  in  an  effort  to  increase 
network  reliability  over  poor  quality  channels,  ARQ  protocols  were  studied  in  the  context  of  coop¬ 
erative  relay  networks  [20]-[22].  In  particular,  [20]  was  among  the  first  such  studies  to  present  a 
general  framework  of  cooperative  ARQ  relay  networks.  It  was  shown  that  cooperative  ARQ  relay 
networks  have  great  advantages  in  terms  of  throughput,  delay,  and  energy  consumption  compared 
to  conventional  multihop  ARQ  networks  in  which  point-to-point  ARQ  links  are  concatenated  to 
form  network  routes.  In  [21],  information-theoretic  analysis  was  developed  and  upper  bounds  for 
the  diversity  order  of  a  decode-and-forward  (DF)  cooperative  ARQ  relay  scheme  were  character¬ 
ized  for  both  slow  and  fast  fading  channels  as  a  means  to  study  the  diversity-multiplexing-delay 
tradeoff.  In  [22],  a  closed-form  expression  of  the  outage  probability  of  the  DF  cooperative  ARQ 
relay  scheme  was  obtained  for  slow  fading  channels,  but,  unfortunately,  the  introduced  approach 
cannot  be  extended  to  fast  fading  channels. 

Outage  probability  is,  arguably,  a  fundamental  performance  metric  for  wireless  ARQ  relay 
schemes  and  so  is  the  diversity  order.  In  this  paper,  we  develop  a  new  analytical  methodology  for 
the  treatment  of  DF  cooperative  ARQ  relay  networks  in  fast  fading  (packet-rate  fading)  channels, 
which  leads,  for  the  first  time,  to  a  closed-form  asymptotically  tight  (as  SNR  ^  oo)  approximation 
of  the  outage  probability.  The  closed-form  expression  shows  that  the  overall  diversity  order  of 
the  DF  cooperative  ARQ  relay  scheme  is  equal  to  2L  —  1,  where  L  is  the  maximum  number  of 
ARQ  retransmissions.  The  achieved  diversity  is  partially  due  to  the  DF  cooperative  relaying  and 
partially  due  to  the  fast  fading  nature  of  the  channels  (temporal  diversity  due  to  retransmissions 
over  independent  fading  channels).  We  note  that  the  diversity  of  the  direct  ARQ  scheme  (without 
relaying)  is  only  L  and  it  is  due  to  the  fast  fading  nature  of  the  channels.  Based  on  the  asymptot¬ 
ically  tight  approximation  of  the  outage  probability,  we  are  able  to  determine  the  optimum  power 
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that  needs  to  be  allocated  at  the  source  and  at  the  relay  of  the  DF  cooperative  ARQ  relay  scheme 
for  any  given  total  transmission  power  budget.  The  optimum  power  allocation  depends  on  the 
variance  values  of  the  channels  involved  and  the  maximum  number  of  retransmissions  allowed  by 
the  protocol.  It  turns  out  that  the  conventional  equal-power  allocation  strategy  is  not  optimum, 
in  general,  and  the  optimum  power  allocation  relies  heavily  on  the  link  quality  of  the  channels 
related  to  the  relay.  Extensive  numerical  and  simulation  results  included  in  this  paper  illustrate 
and  validate  the  theoretical  developments. 

The  paper  is  organized  as  follows.  In  Section  2,  we  describe  briefly  the  DF  cooperative  ARQ 
relay  scheme  and  the  fast  (packet-rate)  fading  channel  model.  In  Section  3,  we  first  develop  two 
useful  lemmas  which  form  the  basis  of  our  analytical  approach.  Then  we  derive  our  asymptotically 
tight  approximation  of  the  outage  probability  of  DF  cooperative  ARQ  relay  scheme.  In  this  section, 
we  also  include  the  outage  probability  expression  of  the  direct  ARQ  transmission  scheme  for  com¬ 
parison  purposes.  In  Section  4,  we  determine  the  optimum  power  allocation  for  the  DF  cooperative 
ARQ  relay  scheme,  and  in  Section  5  we  present  numerical  and  simulation  studies.  Finally,  some 
conclusions  are  drawn  in  Section  6. 

2  System  Model 

We  consider  a  cooperative  ARQ  relay  scheme  with  one  source,  one  relay  and  one  destination  as 
illustrated  in  Fig.  1.  The  DF  cooperative  ARQ  relay  scheme  works  as  follows.  First,  a  data 
packet  of  b  bits  is  encoded  into  a  sequence  of  length  LT,  where  L  is  the  maximum  number  of  ARQ 
retransmission  rounds  allowed  in  the  protocol  and  T  is  the  duration  of  a  single  ARQ  retransmission. 
Then,  the  sequence  comprises  L  different  blocks  each  of  length  T.  In  each  ARQ  retransmission 
round  a  block  of  the  message  is  sent,  so  the  transmission  rate  is  i?  =  h/T.  When  the  source  transmits 
a  block  of  the  message  to  the  destination,  it  is  also  received  by  the  relay.  The  destination  indicates 
success  or  failure  of  receiving  the  message  by  feeding  back  a  single  bit  of  acknowledgement  (ACK) 
or  negtive-acknowledgement  (NACK).  The  feedback  is  assumed  to  be  detected  reliably  both  by  the 
source  and  by  the  relay.  If  an  ACK  is  received  or  the  retransmission  reaches  the  maximum  number 
of  rounds,  the  source  stops  transmitting  the  current  message  and  starts  transmitting  a  new  message. 
If  a  NACK  is  received  and  the  retransmission  has  not  reached  the  maximum  number  of  rounds. 
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Figure  1:  Illustration  of  the  cooperative  ARQ  relay  scheme  with  one  source,  one  relay  and  one 
destination. 


the  source  sends  another  block  of  the  same  message.  If  the  relay  decodes  successfully  before  the 
destination  is  able  to,  the  relay  starts  cooperating  with  the  source  by  transmitting  corresponding 
blocks  of  the  message  to  the  destination  by  using  a  space-time  transmission  scheme  [21]  (e.g.  the 
Alamouti  scheme  [23]).  The  destination  combines  the  received  signal  in  current  round  and  those 
in  previous  rounds  to  jointly  decode  the  data  packet.  After  L  ARQ  retransmission  rounds,  if  the 
destination  still  cannot  decode  the  data  packet,  an  outage  is  declared  which  means  that  the  mutual 
information  of  the  DF  cooperative  ARQ  relay  channel  is  below  the  transmission  rate. 

The  DF  cooperative  ARQ  relay  scheme  can  be  modeled  as  follows.  The  received  signal  yr,m  at 
the  relay  at  the  m-th  (1  <  m  <  L)  ARQ  retransmission  round  can  be  modeled  as 


Vr., 


,m  —  's/ Pshsr^m^s 


(2.1) 


where  Pg  is  the  transmitted  power  of  the  source  signal  Xg,  is  the  coefficient  of  the  source- relay 
channel  at  the  m-th  ARQ  retransmission  round,  and  rjr^m  is  the  additive  noise.  If  the  relay  is  not 
involved  in  forwarding,  the  received  signal  yd,m  at  the  destination  at  the  m-th  ARQ  retransmission 
round  is 

yd,m  ~  \/ Pshgd^m^s  T  'nd,m-i  i^-^) 

where  hgd^m  is  the  source-destination  channel  coefficient  at  the  m-th  ARQ  retransmission  round.  If 
the  relay  receives  the  data  packet  from  the  source  successfully,  it  helps  in  forwarding  the  packet  to 
the  destination  using  the  Alamouti  scheme.  Specifically,  each  block  Xg  of  the  data  sequence  can  be 
considered  as  having  two  parts,  and  Xs,2  (i-e.  Xg  =  [x^q  Xs,2])-  The  relay  forwards  the  block 
in  the  form  of  Xr  =  [— 2  i]  •  received  signal  yd,m  at  the  destination  at  the  m-th  ARQ 
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retransmission  round  can  be  written  as 


yd,m  —  \/ Pshsd,mXs  “1“  \/ Prhrd,mXr  “1“  yd,mi  (^-3) 

where  Pr  is  the  transmitted  power  at  the  relay  and  hrd,m  is  the  channel  coefficient  from  the  relay 
to  the  destination  at  the  m-th  ARQ  retransmission  round.  At  the  destination,  the  message  block 
Xs  can  be  recovered  based  on  the  orthogonal  structure  of  the  Alamouti  code  [21,  23].  The  channel 
coefficients  hsd,m-,  hsr,m  and  hrd,m  are  modeled  as  independent,  zero-mean  complex  Gaussian  ran¬ 
dom  variables  with  variance  and  respectively.  We  consider  a  fast  fading  scenario,  i.e. 

the  channels  remain  fixed  within  one  ARQ  retransmission  round,  but  change  independently  from 
one  round  to  another  (packet-rate  fading).  The  channel  state  information  is  assumed  to  be  known 
at  the  receiver  and  unknown  at  the  transmitter.  The  noise  terms  rjr^m  and  rjd^m  are  modeled  as 
zero- mean  complex  Gaussian  random  variables  with  variance  A/q. 

3  Outage  Probability  Analysis 

3.1  Two  Lemmas 

First  we  develop  two  lemmas  that  will  paly  a  key  role  in  analyzing  the  outage  probability  of  the 
DF  cooperative  ARQ  relay  scheme. 

Lemma  3.1  independent  sealar  random  variables  satisfying  the 

following  properties 

M 

n  <t]=  a-  f{t), 

M 

sl^oo  n  <t]  =  b-  g{t), 

where  di,d2,a  and  b  are  eonstants,  f{f)  and  g{t)  are  monotonieally  inereasing  funetions,  and  f'{t) 
is  integrable,  then 

M  ..t 

<t]=ab-  /  g{x)f'{t-  x)dx.  (3.1) 
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Proof  :  For  any  partition  U  of  the  interval  [0,  t],  i.e.  U  =  {uo,ui,  ...,uj}  with  mq  =  0  and 

uj  =  t,  we  can  obtain  upper  and  lower  bounds  of  the  event  follows: 

J 

i=i 

J 

i=i 

The  upper  and  lower  bounds  in  (3.2)  and  (3.3),  respectively,  are  considered  as  a  union  of  rectangles 
between  uj-i  and  Uj  for  0  <  j  <  J.  First,  let  us  focus  on  the  upper  bound.  The  probability  of  the 
subset  {rij-i  <  Us-i^^,,,^sm  <  ^  {'1’si,...,sm  <  t  ~  %-i}  can  be  calculated  as 

Pr  [uj-i  <  <  Ui,  <t-  %-i] 

=  {Pr  [us,,...,SM  <  ^i]  -  Pr  [Usu...,SM  <  «i-i]}Pr  [r’si,...,SM  <  t  -  uj-i]  .  (3.4) 


Then, 


M 

IT  ■  Pr  K-i  <  risi,...,SM  <  <  t  -  Uj-i]  =  ab-{f{uj)  -  f{uj-i)}g{t  -  uj-i), 

^CXJ 


(3.5) 


and  for  any  monotonically  increasing  functions  f{t)  and  g{t), 

M  J 

sciP[/  n  <t]  <ab-J29{t-  Uj-i){f{uj)  -  f{uj-i)}.  (3.6) 

l<i<Mi=l  j=l 

Similarly, 


M  L 

™fc/  n  ■  Pi'  <t]>ab-^g{t-  Uj){f{uj)  -  (3.7) 

l<i<Mi=l  j=l 

The  above  upper  and  lower  bounds  are  good  for  any  partition  U  =  {uq,  ui, ...,  uj}  over  the  interval 

[0,  f\.  Since  f'{t)  is  integrable,  then,  for  J  ^  oo,  the  sum  terms  in  (3.6)  and  (3.7)  converge  to  the 
same  integral  Jq  g{x)f{t  —  x)dx.  Therefore  we  have  the  result  in  (3.1).  ■ 


We  note  that  the  special  case  of  Lemma  1  with  M  =  1  was  presented  in  [24].  Lemma  3.1  will 
be  used  to  approximate  the  outage  probability  of  the  ARQ  schemes  at  high  SNR  scenario.  In  the 
following,  we  define  a  key  function  Fn{Pi,  ■■■,  (3n',t)  which  will  be  used  to  characterized  the  outage 
probability  of  the  DF  cooperative  ARQ  relay  scheme.  For  any  integer  n  >  2  and  non-zero  constants 
Pi,  P2,  ■■■  ,Pn,  define 

ft  fXn  fX2 

Fn{Pl,...,PnP)=  /  •••/  2^l"i+^2X2+-+/l„a.„d^^dX2---dx„.  (3.8) 

Jo  Jo  Jo 
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The  following  lemma  evaluates  Fn{j3i,  ...,Pn]t)  in  closed-form^. 


Lemma  3.2  For  any  integer  n  >2  and  non-zero  eonstants  /3i,  /32,  •  •  •  ,  fdn,  the  funetion  Fn{j3i, /?„;  t) 
ean  he  ealeulated  as  follows 


6{0,1} 


(3.9) 


where  the  variables  61,82,  G  {0,1},  S  =  {61,62,  ■■■,8n-i\ ,  and  the  eoeffieients  {im,i{S)  ■  1  < 

m  <  n,l  <  I  <  m}  are  speeified  as 


il,liS)  =  i2,2{S)  =  ■■■  =  in,n{S)  =  1, 


and,  for  any  m  =  2,3, ...,  n, 


—  6m— 1  ■  l,z(^))  t  —  1,2,  ...,  m  1. 


Proof  :  We  use  induction  to  prove  the  result  for  any  integer  n  >2.  When  n  =  2,  it  is  easy  to 
see  that 

rt  rx2 


nX2 

2F^^+h2^^dxidx2 

(ln2)-i  f  (ln2)-i  _  (ln2)-i 


/3i  I  /3i  +  P2 

E 


A 


-  1 


(_l)2+5i(in2)-2 


2[E?=i  *2.i(<5)/3i]t  _ 


(3.10) 


(5ie{0,l}  nm=l  [  Z)z=l  *m,z(‘^)A] 
i.e.  the  closed-form  expression  in  (3.9)  is  valid  for  n  =  2.  Next,  we  assume  that  the  result  in  (3.9) 
is  good  for  n  =  k,  where  k  is  any  fixed  integer  greater  or  equal  to  2.  Then,  for  n  =  k  +  1, 

...  2h^-^+h^-^+-+h>^+^-^+idxidx2  •  •  •  dXk+l 

JO 

ft  I-Xk+I  fX2 


n^fc+1 

...  dxidx2  ■■■dxk  2^''+i^'=+i  dxk+i 

Jo 


[  Fk{l3i,...,(3k;  Xfc+i)2^'=+i*''+idxfc+i. 
Jo 


(3.11) 


^Arguably,  the  significance  of  Lemma  3.2  goes  beyond  the  problem  considered  in  this  paper. 
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According  to  the  induction  assumption,  we  have 


r^k+i  f^k  f^2 

FkiPi,...,(3k-,Xk+i)  =  /  •••  / 

Jo  Jo  Jo 

(_l)fc+<5iH  l-(5fc-i  (ln2)“^ 


aS-.O^i  [Er=iW(-5)A] 

6(0,1} 


(2[2:  ik,l  (S)(3i]xk+i  _  _  (3.12) 


Substituting  (3.12)  into  (3.11),  we  obtain 

J^k+l{Pl )  •  •  •)  Pk+1 ) 

(_l)*:+<5iH - ^^fe-i(ln2)“^  /”* 


E 


(5i,...,(5fc_i 

6(0,1} 


n^=l[T.T=l^mAm]  Jo 


^\Yl’i^i'^k,i{^)Pi+Pk+i\xk+i  _  2^'“+i*'=+i^dxfc+i 


_ "I  — \-Sk-i 


(-1) 


A, 5-1  nE.  [E£.  W(<5)A]  ■  [Ef.i4j('S)A  +  A»+ 

6(0,1} 


nE.  [E£.  w(<5)fl]-A+1 

6(0,1} 

^ - \-5k  Qn2)~ 


^2^fc+i 


- 1 


lb.  n^rJ-A[EZl^n.Am] 

6(0,1} 


(3.13) 


where  4  G  {0,1},  ik+i,k+i{S)  =  1  and  ik+i,i{S)  =  6k  ■ik,i{S),l  =  1,2,...,  A:.  Also,  we  have 


k+l 


y^jk+i,i{^)Pi 


1=1 


_  /  Ef=i  ik,i{S)Pi  +  4+1,  if  6k  =  1; 
4+1,  if  4  =  0. 


(3.14) 


Therefore,  the  closed-form  expression  in  (3.9)  is  valid  for  n  =  k  +  1,  which  completes  the  proof. 


3.2  Outage  Probability  of  the  Direct  ARQ  Transmission  Scheme 


For  comparison  purposes,  in  this  subsection,  we  evaluate  the  outage  probability  of  the  direct  ARQ 
transmission  scheme.  The  destination  of  a  direct  ARQ  transmission  scheme  receives  information 
from  the  source  directly,  without  involving  the  relay.  The  mutual  information  between  the  source 
and  the  destination  in  the  m-th  round  of  the  direct  ARQ  transmission  scheme  is 

Isd,m  =  log2  ^1  +  •  (3.15) 

The  total  mutual  information  after  L  ARQ  rounds  is  /{j)*  =  Isd,m-  Thus,  the  outage  proba¬ 

bility  of  the  direct  ARQ  scheme  after  L  ARQ  rounds  is 


pout,L  ^  Pj.  ^jtot  ^  ^ 


(3.16) 


A  closed-form  expression  of  (3.16)  is  not  tractable.  An  approximation  of  the  outage  probability  for 
high-SNR  can  be  obtained  as  a  direct  application  of  Lemma  3.1  with  M  =  1: 


~  ql^R)  ( 


Mo 

(j'i  p 


(3.17) 


where  5'l(-)  is  defined  as 


9n{t)  =  [  gn-i{x)f{t  -  x)dx,  n>l,  (3.18) 

Jo 

with  go{t)  =  1  and  f{t)  =  2*  —  1.  Expressions  (3.17)  and  (3.18)  appeared  also  in  [24].  However,  the 
calculation  of  the  coefficient  gi^R)  that  involves  L  recursive  integrals  was  not  treated  in  [24].  In 
this  paper,  we  develop  a  closed-form  expression  for  the  function  gn{t)  for  any  n  >  1.  Specifically, 
since  gn{t)  =  Jo  9n-i{x)f'{t  —  x)dx  and  f'{t)  =  2dn2,  we  have 

ft  fXn-l  /•X2 


r-t  j'Xn-i  rX2 

9n{t)  =  //  •••/  gi{xi)f{x2  -  Xi)f{x3  -  X2)  ■  ■  ■  f'iXn-1  -  Xn-2)f{t  -  Xn-l)dxidX2- ■ -dXn-l 

JoJo  Jo 

-1  rx2 

■  (2"^i  -  l)(ln2)’^“^2"^2-xi2^3-^2  . . .  2"^"-i-"^"-22*""^"-idxidx2  •  •  •  dx^-i 

Jo 


frr 


ft  fXn-1  fX2 

=  2\ln2)^-^  /  /  ■■■  {!-  2-^i)dxidx2  •  •  •  dx„_i 

JoJo  Jo 


=  2 


(ln2) 


OJO 
n—1  rt 


-  /x”-2(l-2-*+^)dx. 

■  Jo 


{n-2)'.  jQ 
Since  x^~'^  dx  =  and  ([25]) 


n—1 


^71  —  2 


X 


2^  dx  =  -2*  ^ 


(n-2) 


fn—m—l  _|_ 


10  (-ln2)”^(n-m- 1)!  (-In2)"-i^ 

therefore,  a  closed-form  expression  of  gn{t)  can  be  obtained  as  follows 


(n  —  1)! 


m=l 


(n  —  m  —  1)! 


— ^  (m  —  1)! 

-2=1  ^ 

which  can  be  calculated  efficiently. 


^  / _ 1  \n—m 

2‘E  L  .M  (Mn2r-‘+(-ir, 


(3.19) 


3.3  Outage  Probability  of  the  DF  Cooperative  ARQ  Relay  Scheme 


In  this  subsection,  we  derive  the  outage  probability  of  the  DF  cooperative  ARQ  relay  scheme  under 
packet-rate  fading  conditions.  In  the  DF  cooperative  ARQ  relay  scheme,  if  the  relay  decodes  the 
message  from  the  source  correctly,  say,  at  the  k-th  round,  then  at  the  {k  +  l)-th  round,  the  relay 
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starts  forwarding  appropriate  ARQ  blocks  to  the  destination.  Let  {Tr  =  k}  denote  the  event  of 
successful  message  decoding  by  the  relay  at  the  fc-th  round  and  subsequent  ARQ  block  forwarding 
at  the  {k  +  l)-th  round.  Let  denote  the  conditional  probability  that  the  destination  decodes 

the  message  unsuccessfully  after  L  ARQ  retransmission  rounds  given  that  the  event  {R.  =  k} 
occurred.  In  other  words,  P^Lj.  denote  the  outage  probability  at  the  destination  despite  the  fact 
that  the  relay  started  forwarding  at  the  {k  +  l)-th  round.  Therefore,  the  outage  probability  of  the 
DF  cooperative  ARQ  relay  scheme  after  L  ARQ  retransmission  rounds  can  be  written  as 

L 

pout,L  ^  ^  ^  (3  20) 

k=l 

First,  we  calculate  the  probability  of  the  event  {R.  =  /c}  ,  i.e.  Pr  [T^  =  k\.  The  mutual  infor¬ 
mation  between  the  source  and  the  relay  in  the  m-th  ARQ  round  is 

Isr,m  =  log2  +  ~J^\hsr,m\"^  ■  (3-21) 


We  note  that  the  channels  change  independently  over  each  ARQ  retransmission  round  in  a  fast 
fading  scenario,  so  the  mutual  information  of  fading  channels  can  be  viewed  as  a  sum  of  independent 
random  variables.  The  probability  that  the  relay  decodes  the  message  successfully  at  the  first  round 
{Tr  =  1)  is 

Pr  [Tr  =  1]  =  Pr  [hr,i  >R]=  exp  ■  (3.22) 

V  ^sr  / 

For  any  Tr  =  k,  k  =  2,  3, ...,  L  —  1,  we  have 


Pr  [Tr  =  k]  =  Pr 


fc-i 


^  ^  Isr,m  ^  Pi  ^  ^  ksr,m  P  P 


Lm=l 

fc-1 


=  Pr 


^  ^  Isr,m  ^  P 


m=l 


gk-i{P) 


Mo 

a'lrPs 


m=l 

—  Pr 

k-l 


^  ^  ksr,m  ^  P 


m=l 


-  9k{P) 


Mo 

C^srPs 


(3.23) 


where  gk-i{-)  and  gk{-)  are  specified  in  (3.19)  in  the  previous  subsection.  The  approximation  in 
(3.23)  is  obtained  by  applying  Lemma  1  with  M  =  1.  Finally,  if  Tr  =  L,  we  have 


Pr  [Tr  =  L]=Pr 


L-l 

E- 

|_m=l 


<  p 


'gL-i{P) 


Mq 

(^IrPs 


L-l 


(3.24) 
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To  summarize,  the  probability  Pr  [T,.  =  k\  can  be  given  by  the  following  branch  function 


2^-1  .  M)  \ 
Csr  Ps  J  ’ 


k  =  1; 


Pr(r,  =  fc|~l  gt-iiR)  {^y  '-gk{R){^y\  2<k<L-l; 
gL-,(S)(£fry~‘,  K  =  L. 


(3.25) 


Next,  we  calculate  the  conditional  outage  probability  when  the  relay  decodes  correctly 

at  the  fe-th  round  and  starts  forwarding  at  the  {k  +  l)-th  round.  This  is  done  by  the  following 
theorem. 


Theorem  3.1  The  conditional  outage  probability  1  <k  <  L,  is  given  by 


bk(R)  (  pQ  \^  (  pQ  ^  -x  ^  u  ^  T  1. 


(3.26) 


gm(A, 


k  =  L, 


where  gk{-)  are  specified  in  (3.19),  and 


bk{t)=  /  gk{x)qL_k{t- x)dx, 
Jo 


(3.27) 


in  which  qi{t)  =  (2*  —  1)^  and  for  any  2  <  n  <  L  —  1,  qn{t)  is  given  by 


Qn{t)  =  {-2ln2Y  (-l)"i+'''+""-i2(^+""-i)*{F„_i(l-ai,ai-a2,---  ,an-2-an-i;t) 


ai,---  iQn-l 

6(0,1} 


2‘Fn—li,  CXi,CXi  02)  ■  ■  ■  1  Oin—2  Oin—lit) 

+  Pn-i(— l  —  «i,  01  —  02)  ••  •  )  an_2  — On-i;  t)}.  (3.28) 


The  function  Fn-i{  •  ]t)  is  specified  in  Lemma  3.2. 


Proof  :  When  the  relay  cooperates  with  the  source  by  jointly  sending  a  message  block  via  the 
Alamouti  scheme,  the  mutual  information  of  the  cooperative  channels  in  the  m-th  ARQ  round  is 


given  by  [21] 


f  1 2  1 2^ 

Isrd,m  —  log2  (  ^  T  1 


(3.29) 


Thus,  with  L  ARQ  rounds,  the  total  mutual  information  is 


jtot  _ 

^d,Tr=k  — 


'Ylm=l  ^sd,m  +  X]m=fc+1  ^srd,mi  1  ^  k  <  L] 


(3.30) 


k  =  L. 
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The  above  mutual  information  is  based  on  the  assumption  that  the  channels  change  independently 
over  each  ARQ  retransmission  round,  so  the  mutual  information  of  fading  channels  can  be  viewed 
as  a  sum  of  independent  random  variables.  We  also  note  that  when  =  L,  the  relay  has  no  chance 
to  cooperate  since  the  source  starts  sending  a  new  packet. 

The  conditional  outage  probability  P^l!'Lk  be  evaluated  as 


Pt=k  =  Pr  <  R]  ■ 


(3.31) 


When  Tr  =  L,  the  conditional  outage  probability  is  reduced  to  the  direct  ARQ  scenario  and  it  is 


given  by 


P^:U  =  Pt 


(3.32) 


In  the  following,  we  calculate  the  conditional  outage  probability  (3.31)  for  any  Tr  =  k,  k  = 


1, 2, ...,  L  —  1.  For  simplicity  in  presentation,  we  introduce  the  following  notation.  Let 

|'log2  (l  +  ,  1  <  m  <  A:; 


[log2  (l  +  Sl\hsd,7n\‘^  +  S2\hrd,7n\‘^)  ,  k  +  1  <  m  <  L, 

where  si  =  Pg/Mo  and  S2  =  Pt/Nq.  Then,  the  total  mutual  information  can  be  written  as 

L 

^d°Tr=k  =  ^ 


(3.33) 


(3.34) 


We  note  that  for  any  1  <  m  <  A:,  \hsd,m\‘^  is  an  exponential  random  variable  with  parameter  cr^^. 


Thus, 


r  2*  —  l" 

lim  Si  •  Pr  [um  <  t]  =  lim  si  •  Pr  \hsd,raf  <  - 

51^00  51^00  5]^ 


(2‘-l),  m  =  l,2---,A.  (3.35) 


Since  Um,  1  <  m  <  A:,  are  independent  random  variables,  by  applying  Lemma  1  with  M  =  1 


recursively,  we  have 


lim  Si  •  Pr  ^  Um  <  A  =  (  — )  9k{t) 


(3.36) 


where  gk{t)  is  given  in  (3.19).  For  any  m,  k  +  1  <  tu  <  L,  Um  involves  the  sum  of  two  independent 
exponential  random  variables  |  hsd,m  \  ^  and  |  hrd,m  \  ^  with  parameters  and  ,  respectively,  and 
the  distribution  of  Um  can  be  specified  as 


Pr  [Um  <  A]  =  Pr  [si\hsd,ra?  +  S2\hrd,ra?  <  2*  -  l] 


1-1 


_  ^2_ 


1 - ^  2  exp 


-exp  (-^2^),  if;^^;^, 

1  \  ^rd  2  /  '^rd  ‘^sd 


(3.37) 
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for  any  m  =  /c  +  1, L.  Thus,  for  any  m  =  k  +  1, L,  we  have 


•  Pr  [u^  <t\  =  ^  o  _2  -  1)^- 


(3.38) 


Let  qQ{t)  =  1  and  p{t)  =  (2*  —  1)^,  then  p' {t)  =  2(2^*  —  2*)ln2.  Since  Um,  A:  +  1  <  m  <  L,  are 
independent  to  each  other,  by  applying  Lemma  3.1  with  M  =  2  recursively,  we  can  show  that  for 
any  n  =  1, 2, ...,  L  —  k, 


Ji^  (S1S2)"'  •  Pr  ^  Um  <t 


.m=/c+l 


in  which 


qn{t)=  qn-i{x)p'{t  -  x)dx,  n  =  1,2,  ...,L  -  k. 

Jo 


(3.39) 


(3.40) 


We  can  see  that  when  n  =  1,  qi{t)  =  {2^  —  1)'^ .  But  for  larger  n,  the  calculation  of  qn{t)  is  involved. 
Based  on  Lemma  3.2,  a  closed-form  expression  for  qn{t)  can  be  obtained  as  follows: 

rt  f-Xn-l  l'X2 

qn{t)  =  //  •••/  qi{xi)p  {X2  -  Xl)p  {X2,  -  X2)  ■  ■  ■  p  {Xn-l  -  Xn-2) 

JoJo  Jo 


xp'{t  -  Xn-l)dxidX2  ■  ■  ■  dXn-l 


ft  fXn-l  fX2 

/  /  •••/  (2^1  -  l)2(21n2)’^-i  TT 

JoJo  Jo  ^  J 


x(2*  3;„-i  _  *idxidx2  •  •  -dXn-l 


(-21n2)'"“^  E 


G{0,1} 

^  r^n—i  r^2 


rt  pXn-l  rX2 

X  •  •  7  (2*^  -  1)^  •  IT  2("™-i““’")*™dxidx2  •  •  •  dxn-i 

JoJo  Jo  ^^2 

(-21n2)"-i  (-l)“i+-+“-i2P+“"-P7Fn_i(l-ai,/32,...,/3n-i;t) 


6(0,1} 


2Tn_l  (  Cti ,  (52 1  •  •  •)  Pn—lj  a)  T  Pn— 1  (  1  Cli ,  P2  j  ■■■j  Pn—lj  A)  }  > 


(3.41) 


where  ^2  =  0^1  —  02;  f5o  =  ot2  —  03,  ...,  f5n-i  =  Q;n-2  —  cin-i,  and  P’„_i(  •  ;  t)  is  defined  in  (3.9)  (a 
closed-form  expression  for  Fn-i{  •  ;  A)  is  given  in  Lemma  3.2).  We  note  that  in  the  case  where  f3i  is 
zero  for  some  i,  the  non-zero  condition  in  Lemma  3.2  is  not  satisfied.  In  such  case,  we  may  evaluate 
F„_i(  •  ;  t)  by  applying  Lemma  3.2  with  (5i  =  Si  where  Si  is  sufficiently  small  (i.e.  £i  0)  (Since 
the  function  Fn{(5i,  ...,(5n',t)  defined  in  (3.9)  is  continuous  in  terms  of  /3i,  2  <  i  <  n  —  1,  so  is  the 
closed- form  expression  in  Lemma  3.2). 
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According  to  the  result  in  (3.39)  with  n  =  L  —  k,  we  have 


lim  (siS2)^  ^  •  Pr 

5^— >00  ^  ' 

l<i<2 


Urn  <t 

\jm=k+l 


L-k 


QL-kit), 


(3.42) 


where  can  be  calculated  specifically  based  on  (3.41).  Combining  (3.36)  and  (3.42),  and 

applying  Lemma  1,  we  obtain 


lim  sf  Sn  ^  •  Pr 

Si— >oo  ^  ^ 

l<i<2 


k  L 

^  ^  ksd,m  ^  ^  Isrd,m  ^  R 
m=k+l 


m=l 


=  bk{R)  ( 


a 


sd 


^^rd 


L-k 


(3.43) 


where 


hit)  =  [  gk{x)qL_k{t  -  x)dx.  (3.44) 

Jo 

Since  si  =  Pg/J^o  and  =  Pr/A/o,  so  for  any  Tr  =  k,  k  =  1,2, ...,  L  —  1,  the  conditional  probability 
(3.31)  can  be  asymptotically  approximated  as 


r  k 


Pn=k  =  Pr 


^  ^  ksd,m  P  ^  ^  J-srd,m  P 


R 


m=l  m=k+l 

which  completes  the  proof  of  the  theorem. 


bk(R)  f  Mo  V  r  M 


^IdPr 


L-k 


(3.45) 


Finally,  based  on  the  probability  Pr  \Tr  =  k]  in  (3.25)  and  the  conditional  outage  probability 
PjJ^k  ™  Theorem  3.1,  we  can  obtain  the  outage  probability  for  the  DF  cooperative  ARQ  relay 
scheme  as  follows 


P' 


>out,L 


L-1 

E 


him 

2L-k 


gk-i{R)  -  gk{R) 


k=l 

+  gL{R)gL-i{R) 


A/q 

fj2  p 


A/q 

Cr^rPsA 

A/i) 

a'l^Ps 


A/q 

L-1 


A/q 


L-k 


A/q 

MsrPs 


k-l 


(3.46) 


Furthermore,  we  note  that  the  term  gk{R)  would  be  much  smaller  than  gk-i{R)  at  high  SNR 

<7  grP^s 

so  the  asymptotic  outage  probability  in  (3.46)  can  be  further  simplified  as 


,onLL  ^bk{R)gk-i{R)  (  A/))  \V  A/),  A/))  V"' 

y^'dPs)  WrdPr)  VlPg) 


(3.47) 


where  biiR)  =  gL{R)-  Simulation  studies  presented  later  in  this  paper  will  illustrate  the  tightness 
of  the  outage  probability  at  high  SNR. 

Based  on  the  above  asymptotic  outage  probability,  we  observe  that  the  term  (  )  in  (3.47) 

contributes  a  diversity  order  L  in  the  asymptotic  outage  performance,  which  is  due  to  the  fast 
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fading  nature  of  the  channels.  The  term  (  ^  °  ^  p  )  contributes  an  overall  diversity 

order  (L  —  k)  +  {k  —  1)  =  L  —  1  which  is  due  to  the  cooperative  relaying.  Thus,  the  asymptotic 
outage  probability  of  the  DF  cooperative  ARQ  relay  scheme  has  an  overall  diversity  order  2L  —  1. 
In  the  case  of  the  equal  power  allocation,  i.e.  Pg  =  Pr  =  P,  the  contribution  of  the  diversity  order 
in  the  outage  probability  is  more  evident.  We  recall  that  the  diversity  order  of  the  direct  ARQ 
transmission  scheme  is  only  L,  which  is  due  to  the  fast  fading  nature  of  the  channels,  and  it  is 
much  smaller  than  that  of  the  DF  cooperative  ARQ  relay  scheme. 


4  Optimum  Power  Allocation  for  the  DF  Cooperative  ARQ  Relay 
Scheme 


In  this  section,  we  derive  the  asymptotically  optimum  power  allocation  strategy  for  the  DF  cooper¬ 
ative  ARQ  relay  scheme  based  on  the  approximation  of  the  outage  probability  that  was  presented 
in  the  previous  section.  Without  loss  of  generality,  we  denote  the  total  transmission  power  as 
Ps  +  Pr  =  2P,  where  Pg  and  Pr  are  the  power  used  by  the  source  and  the  relay,  respectively.  Then, 
for  any  given  total  transmission  power  2P,  we  try  to  determine  optimum  power  Pg  and  Pr  in  order 
to  minimize  the  asymptotic  outage  probability. 


Let  A  denote  the  ratio  of  the  source  power  Pg  to  the  total  transmission  power,  i.e.  A  =  ^. 
Then  0  <  A  <  1  and  Pr  =  (1  —  X)2P.  The  asymptotic  outage  probability  of  the  DF  cooperative 
ARQ  relay  scheme  can  be  written  as 


pout^L 


A/i) 


2L-1  L 

Y,hk{R)gk-i{R) 

k=l 


1 


(4.1) 


We  try  to  find  the  optimum  power  ratio  A  (0  <  A  <  1)  such  that  the  asymptotic  outage  probability 


A 


is  minimized.  Let  Ak{R)  =  bk{R)gk-i{R) 


G{\)  = 


and 

Ak{R) 


^  A^+^-i(l  -  A) 
Then  the  optimization  problem  can  be  formulated  as  follows 


L-k  • 


(4.2) 


m^n  G'(A) 

s.t.  0  <  A  <  1.  (4.3) 
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(4.4) 


The  optimum  value  of  A  (i.e.  \opt)  satisfies  the  following  equation 

dG{\)  wpj  -{L  +  k-l)  L-k 


or  equivalently 

|-(L  +  t  -  1)  +  (i  -  k)  ( V)'  '}  =  "■ 


(4.5) 


Equation  (4.5)  can  be  easily  solved  by  the  Newton  method.  Looking  closely  at  equation  (4.5),  we 
observe  the  followings.  Since  Ak{R)  is  positive,  then  must  be  less  than  1,  otherwise  the  left- 
hand  side  of  (4.5)  is  negative.  It  is  thus  implied  that  A  >  ^,  i.e.  Pg  >  P  and  Pr  <  P  which  means 
that  we  should  allocate  more  power  at  the  source  and  less  power  at  the  relay.  It  also  shows  that  the 
equal  power  allocation  scheme  that  assigns  equal  power  to  the  source  and  the  relay  is  not  optimum 
in  general.  On  the  other  hand,  for  any  given  transmission  rate  R,  the  parameters  Af^^R)  in  (4.5) 
depend  only  on  and  which  are  the  variance  values  of  the  source-relay  and  relay-destination 
channel  links,  respectively.  Thus,  the  asymptotic  optimum  power  ratio  Xopt  depends  only  on  the 
the  variance  of  the  source- relay  and  relay-destination  channels,  and  not  on  the  source-destination 
channel  link.  A  similar  observation  was  reported  in  [11]  where  the  optimum  power  allocation 
between  the  source  and  the  relay  was  determined  based  on  the  analysis  of  the  symbol-error-rate 
performance. 


As  an  example,  we  study  the  case  where  L  =  2.  In  this  case,  equation  (4.5)  is  reduced  to 


3A2iR) 

so  the  optimum  power  ratio  is 


1  -  A 


+  2Ai(i?) 


1  -  A 


-  Ai(i?)  =0, 


(4.6) 


A  = 


1  -|-  -i  / 1  -|-  3 


A2{R) 

AiiR) 


2  -|-  -i /l  -|-  3 


A2{R) 

AiiR) 


(4.7) 


Thus  the  corresponding  optimum  power  allocation  at  the  source  and  at  the  relay  is  given  by 


Ps  = 


1  -|-  -v  / 1  -|-  3 


A2iR) 

AiiR) 


Pr  = 


p, 


p. 


(4.8) 

(4.9) 


In  (4.8)  and  (4.9),  ^  Thus,  the  optimum  power  allocation  dependents  only 

on  the  ratio  of  the  variance  of  the  source-relay  channel  and  the  variance  of  the  relay-destination 
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channel,  which  is  consistent  with  our  observation  that  we  based  on  equation  (4.5).  Also,  from  (4.8) 
and  (4.9),  we  can  see  that  P  <  Pg  <  2P  and  0  <  Pr  <  P,  i.e.  we  should  allocate  more  power 
at  the  source  and  less  power  at  the  relay  to  optimize  the  overall  performance  at  the  destination. 
Furthermore,  (4.8)  and  (4.9)  imply  that  if  the  relay  is  located  close  to  the  source,  i.e.  ^ 

Pg  goes  to  |P  and  Pr  goes  to  |P.  On  the  contrary,  if  the  relay  is  located  close  to  the  destination, 
i.e.  a^r  ^  <7^^,  Pg  goes  to  2P  and  Pr  goes  to  0,  which  means  that,  in  this  case  we  should  allocate 
most  of  the  power  (2P)  at  the  source.  The  latter  is  is  reasonable  since  the  cooperative  role  of  the 
relay  is  minor  in  this  case. 

5  Simulation  Results 

In  this  section,  we  present  numerical  and  simulation  studies  that  compare  the  performance  of  the 
DF  cooperative  ARQ  scheme  with  that  of  the  direct  ARQ  scheme.  In  all  studies,  the  variance  of 
the  channel  hij{{i,j)  G  {s,d),{s,r),{r,d)}  is  assumed  to  be  afj  =  (i~T,  where  dij  is  the  distance 
between  two  nodes  and  ^  is  the  path  loss  exponent  which  is  assumed  to  be  /x  =  3  in  a  typical 
fading  environment.  We  assume  that  the  source-destination  distance  is  dgd  =  10  m  and  the  relay  is 
located  in  the  midpoint  between  the  source  and  the  destination.  We  consider  a  target  transmission 
rate  of  i?  =  2  bits/s/Hz. 

Figs.  2,  3  and  4  illustrate  the  performance  curves  when  the  maximum  number  of  ARQ  retrans¬ 
mission  rounds  is  L  =  2,  3  and  4,  respectively.  In  these  figures,  equal  power  allocation  is  assumed, 
i.e.  Pg  =  Pr  =  P.  All  figures  show  that  the  proposed  theoretical  approximation  of  the  outage 
probability  of  the  DF  cooperative  ARQ  relay  scheme  is  tight  for  high  SNR  values  while  it  is  less 
tight  for  low  SNR  values.  For  example,  in  Fig.  2  (L  =  2),  the  analytical  approximation  curve  is 
almost  indistinguishable  from  the  simulated  curve  for  all  outage  performance  levels  below  10“^. 
Moreover,  the  larger  the  number  of  ARQ  retransmission  rounds,  the  higher  the  diversity  order  of 
the  DF  cooperative  ARQ  relay  scheme.  This  observation  is  consistent  with  the  theoretical  result 
that  the  diversity  order  of  the  DF  cooperative  ARQ  relay  scheme  increases  as  the  number  of  ARQ 
retransmission  rounds  increases. 

All  Figs.  2-4  show  that  the  DF  cooperative  ARQ  relay  scheme  significantly  outperforms  the 
direct  ARQ  scheme.  At  an  outage  performance  of  10“^,  the  performance  of  the  DF  cooperative 
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ARQ  relay  scheme  is  about  8dB  better  than  that  of  the  direct  ARQ  scheme.  For  the  same  maxi¬ 
mum  number  of  retransmission  rounds  L,  the  DF  cooperative  ARQ  relay  scheme  exhibits  a  higher 
diversity  order  than  the  direct  ARQ  scheme  which  is  consistent  with  our  theoretical  developments, 
i.e.  the  DF  cooperative  ARQ  relay  scheme  has  diversity  order  2L  —  1,  while  the  direct  ARQ  scheme 
has  diversity  order  only  L. 

In  Fig.  5,  we  plot  the  power  allocation  optimization  function  G(A)  in  terms  of  A  (0  <  A  <  1)  for 
the  cases  where  L  =  2,  3  and  4,  respectively.  We  assume  that  the  quality  of  the  source-relay  link  is 
the  same  as  that  of  the  relay-destination  link,  i.e.  We  observe  that  the  optimum  power 

ratio  A  is  about  0.8  for  the  three  cases.  More  precisely,  from  the  numerical  results,  the  optimum 
power  ratios  are  A  =  0.8203,  A  =  0.7969  and  A  =  0.7838  for  L  =  2,  3  and  4,  respectively.  It  appears 
that  the  optimum  power  ratio  decreases  gradually  when  the  maximum  number  of  retransmission 
rounds  is  increased.  Furthermore,  the  optimum  power  ratio  is  much  larger  than  1/2  which  is 
consistent  with  our  analysis  that  we  should  allocate  more  power  at  the  source  and  less  at  the  relay. 

In  Figs.  6-8,  we  compare  the  performance  of  the  DF  cooperative  ARQ  relay  scheme  with  op¬ 
timum  power  allocation  and  with  equal  power  allocation  for  L  =  2,  3  and  4,  respectively.  Both 
simulation  and  numerical  approximation  curves  are  included.  In  Fig.  6  (L  =  2),  the  optimum 
power  allocated  at  the  source  and  the  relay  is  set  to  the  numerical  values  obtained  from  Fig.  5,  i.e. 
PsI^P  =  0.8203  and  Pr/2P  =  0.1797,  respectively.  We  observe  that  the  performance  of  the  DF 
cooperative  ARQ  relay  scheme  with  the  optimum  power  allocation  is  about  1.5dB  better  than  that 
of  the  scheme  with  the  equal  power  allocation.  In  Fig.  7  (L  =  3),  the  optimum  power  allocation  is 
Ps/2P  =  0.7969  and  Prj2P  =  0.2031,  and  the  corresponding  curve  shows  a  performance  improve¬ 
ment  of  1.25dB  compared  to  the  equal  power  allocation  case.  In  Fig.  8  {L  =  4),  the  optimum 
power  allocation  is  Ps/2P  =  0.7838  and  Pr/2P  =  0.2162.  We  see  that  the  performance  curve  of 
the  DF  cooperative  ARQ  relay  scheme  with  optimum  power  allocation  also  exhibits  performance 
gains  of  about  1.25dB  compared  to  the  performance  of  the  equal  power  allocation  scheme. 

6  Conclusion 

In  this  paper,  we  developed,  for  the  first  time,  a  closed-form  asymptotically  tight  (as  SNR  ^  oo) 
approximation  of  the  outage  probability  for  the  DF  cooperative  ARQ  relay  scheme  under  fast 
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fading  conditions.  The  closed- from  expression  provides  significant  insight  into  the  merits  of  DF 
cooperative  ARQ  relaying  relative  to  the  direct  ARQ  scheme  in  fast  fading  scenarios  and  shows 
that  the  cooperative  scheme  achieves  diversity  order  equal  to  2L  —  1  while  the  diversity  order  of 
the  direct  scheme  is  only  L.  Simulation  and  numerical  studies  illustrated  that  the  closed-form 
approximation  of  the  outage  probability  is  tight  at  high  SNR.  Based  on  the  asymptotically  tight 
approximation  of  the  outage  probability,  we  were  able  to  determine  the  optimum  power  allocation 
strategy  for  the  DF  cooperative  ARQ  relay  scheme.  It  turns  out  that  equal  power  allocation  is  not 
optimum  in  general  and  that  the  optimum  power  allocation  strategy  depends  on  the  link  quality  of 
the  channels  related  to  the  relay.  It  is  also  clear  that  we  should  allocate  more  power  at  the  source 
and  less  at  the  relay.  Further  numerical  and  simulation  studies  illustrated  the  performance  gains  of 
the  DF  cooperative  ARQ  relay  scheme  with  optimum  power  allocation  relative  to  the  equal  power 
allocation  scheme. 
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P/Wo(dB) 

Figure  2:  Outage  probability  of  the  direct  and  DF  cooperative  ARQ  schemes  (L=2). 


P/Ai)(dB) 

Figure  3:  Outage  probability  of  the  direct  and  DF  cooperative  ARQ  schemes  (L=3). 
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Figure  4:  Outage  probability  of  the  direct  and  DF  cooperative  ARQ  schemes  (L=4). 


Figure  5:  Optimum  power  ratio  A  for  the  DF  cooperative  ARQ  scheme.  When  L=2,  3  and 
4,  the  asymptotic  power  allocation  is  A  =  0.8203,  A  =  0.7969  and  A  =  0.7838,  respectively. 
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Figure  6:  Outage  probability  of  the  DF  cooperative  ARQ  scheme  with  equal  and  optimum 
power  allocations  (L=2). 


P/NoidB) 

Figure  7:  Outage  probability  of  the  DF  cooperative  ARQ  scheme  with  equal  and  optimum 
power  allocations  (L=3). 
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P/Wo(dB) 

Figure  8:  Outage  probability  of  the  DF  cooperative  ARQ  scheme  with  equal  and  optimum 
power  allocations  (L=4). 
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